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Abstract 

In frame of a phenomenological approach based on Compton-like Feynman amplitudes, we study 
multi-pion production in antiproton nucleon collisions. The main interest of this reaction is related 
to the possibility to study the properties of the presumable NN atom and of its resonances. For 
the case of formation of a scalar or pseudoscalar resonant state, with I G (J P ) = 1~(0 =1: ), + (0~) 
numerical results are obtained. The differential cross section in an experimental set-up where the 
pions invariant mass is measured, is given with explicit dependence on the lepton pair and pions 
invariant mass. 
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In a previous paper [lj we calculated the differential cross section for single pion produc- 
tion in the reactions p + N — > tc + £ + + £~, and showed that such reactions are measurable at 
upcoming facilities, bringing unique information on electromagnetic and axial nucleon form 
factors in the unphysical region of time-like momentum transfer squared, q 2 . 

The purpose of the present work is to study the resonant production of a number of 
pions accompanied by a lepton pair, in the antiproton nucleon annihilation process. We limit 
our discussion to quantum numbers of the excited resonances: I G (J P ) = l + (0 ± ) and + (0 ± ). 
The existence a NN bound state and of its resonances was predicted in a series of papers, 
see Ref. . A partial wave analysis performed in [3( in the range 1960 < M* < 2410 

combines data in different decay channels 3n°, tt t], 7r°r]r]' and reveals a series of resonant 
states in the pp system. Some of the states were precisely identified, in particular two 
states with quantum numbers I(J PC ) = 1(0 _+ ), with masses M* = 2360 and 2070 MeV, 
respectively and widths about 300 MeV. Higher spin states were also identified. Such states 
show similar masses as dibaryon states, except that they exhibit large widths, which are 
explained by their decay through hadronic states, which is impossible for dibaryons 4j. 

NN state has hydrogen-like atomic structure. Similarly to positronium, in addition 
to strong interaction, QED effects are present, due to the large masses (small distances) 
involved. Special attention is payed to resonances with masses ~ 2M, where M is the 
nucleon mass. The resonances with zero orbital momentum are expected to have large 
width, whereas higher values of orbital momentum lead to a smaller width given by [2|: 

100 

T ~ MeV 1) 

~ [(2£ + l)!!] 2 ' 1 ' 

which gives, for example, V = 100 MeV for £ = 0, T = 10 MeV for £ = 1, and T = 0.5 MeV 
for £ = 2. 

Experimental evidence of such series of dibaryon resonances exist {4], although it is still 
controversial. An investigation of this problem will be possible at the upcoming FAIR facility 
[5I, where the planned antiproton beams render possible the measurement of the reactions 
investigated in the present paper. 

The emission of a lepton pair permits to select the appropriate kinematics adapted to 
the excitation of such resonances, when the total pN center of mass (CMS) energy exceeds 
the mass of the resonances, due to a known mechanism, called 'return to resonances' jf|: 
the lepton pair carries away the extra energy and momentum, providing the condition of 
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FIG. 1: Feynman diagrams for the reaction + n^-vr: lepton pair emission form the 

antiproton (a) and from the proton (b). 

exciting the resonances. Experimentally it manifests through a deformation of the Breit 
Wigner distribution, yielding a 'radiative' tail. For the excitation of a narrow resonance, the 
emission of soft quanta alters the shape of the resonance curve: the height of the resonance 
curve decreases and a radiative tail arises. Indeed, for total energy W, higher than the 
resonance mass M' , the most favourable mechanism is when the energy " excess" W — M' is 
absorbed by the photon emission (with energy u = W — M'). Here the small extra factor a 
is compensated by the resonance denominator. 

Let us consider the annihilation of an antiproton and a proton in a lepton pair, and a 
resonant state which subsequently decays into a set of n n pions: p(pi) + pfa) —> 7*(<z) + 
R(J P , M') — > i + (p + ) + £~{p-) + tti(<?i) + ••• + 7T nx (g„ w ). Such state can have the quantum 
number of scalar, pseudoscalar, axial or vector particle. The Feynman diagrams for such 
reactions are plotted in Fig. [TJ We neglect direct photon emission from the resonant state, 
assuming that it is (indirectly) included via the width of the resonance. 

To describe the kinematics of the process it is convenient to introduce two kinematical 
variables: the invariant mass squared of the lepton pair 



Q 2 = (P+ + P-) 2 



(2) 



and the one of the pions: 




(3) 
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One can write the phase volume as: 

d 3 p + d 3 p- d 3 q 



with 



/*-. = / ( Pl + n - P+ -v- - E,j ^^nt.- w 

= J Jr, J dr%-> J dT Q , (4) 

y rfr Qg = y d A qd A Q5\ Pl +P2 - q -Q). 

— * 

It is convenient to calculate this last quantity in CMS, where <f = — Q, pi = — p*2- In this 
reference frame we have: 

j d A q = l -dq \q\dq 2 dO q = J dT Qq , (5) 

where dO q = 2nd cos 9 q is the element of solid angle, integrated on the azymuthal angel (f>. 
Then, using the relations y/s = q + Q ; (s 1 — q 2 )/ 'y/s = Qo — go? an d s — (pi + P2) 2 , one 
finds: 

f 7T f^-V^ dSi 

/ dT Qq = - / dq 2 dcos9 q y/A(s,s 1 ,q 2 ), A(a,b,c) = a 2 + b 2 + c 2 -2(ab + ac + bc), 

J ^ J n ^ m ^ S 

(6) 

with 9 q = q,p\. The function A(s, s±, q 2 ) must be positive in the physical region. 

Let us restrict our considerations to scalar J p = + and pseudoscalar J p = 0~ states. 
The corresponding matrix elements are: 

Mies = ^V p (q 2 )J,v( Pl )0;u( P2 ) - _ — (7) 

where M' and V are the mass and the width of the resonance, V p (q 2 ) = m 2 /(m 2 — q 2 ) 
corresponds to the coupling of the photon to the p meson (m p is the p meson mass), when 
interacting with a nucleon, J^{q) = v(p + )j^u(p^) is the leptonic current and O 1 ^ 
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(where 7W = j 5 for a pseudoscalar resonant state and 7W = 1 for a scalar state) contains 
the hadronic structure: 

= m 2 )l, + ^-(^ - 7m9), (9) 

which is parametrized in terms of the Pauli and Dirac proton form factors, F\ and F2 (in 
time-like region). 

Let us note that the hadron current = v(pi)0^u(p2) obeys the gauge invariance: 

Ji% = 0. 

In Eq. ([7]), the coupling constant g% is related to the width of the decay channel: 

^ = ^gljdv^\ (io) 

where T' n ^ is the width of the decay of the pp state to pions, g\ represents the coupling 
constant at the pp vertex to the resonant state. When the resonance is heavier than twice 
the nucleon mass its value can be related to the pp decay width: 

/ i\2 _ 87tT pp (11) 
{9l) ~ M'Ai ' [ ' 



with At = f3 for J p = 0", and A, = (3 3 for J p = 0+ and (3 = a/1 - (4M 2 /M /2 ). 
The corresponding differential cross sections are: 

da' _2a 2 T' n r pp y/A( s ,s u q*) 1 



dq 2 d Sl 3vr (si - M' 2 ) 2 + M' 2 Y 2 s 2 r q 2 

and 



V u for M' > 2M, (12) 



da* a 2 M'T' n ( 9i ) 2 JA(s, Sl ,q 2 ) . , . 

V -T>i, for M < 2M, (13) 



dq 2 ds 1 12ttV ( Sl - M' 2 ) 2 + M' 2 T 2 s 2 r 



with r = a/1 - (4M 2 /s) and 

2* = / dcos0 fl (^-^S», (14) 



and 



S$ = -Tr^ - M)0;(p 2 + M)Ol*. (15) 



Two dimensionless variables, X = q 2 /s and Y — si/s can been defined, using the invariant 
mass of the lepton pair, q 2 , Eq. (j2J), and the invariant mass of the pions, Si, Eq. ([3]). The 
kinematical region scanned by the reaction, is defined by the conditions A(s, sx, q 2 ) > and 
it corresponds to Y < (1 — \fX) 2 , as it is shown in Fig. [2j 
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FIG. 2: Allowed kinematical region, as a function of X = q 2 /s and Y = si/s (shaded area) 



The angular integration entering in D l , Eq. (JHJ), can be calculated with the help of the 
following relations: 



d cos 9 q 

d COS q 

di 



d cos 6 q 

d cos 
do 



2s 



q 2 s Sl + M 2 A{s,s 1 ,q 2 )' 



r^A(s,S!,q 2 ) 



log 



q 2 + Si - s + r v /A(s, Si, g 2 



<7 2 + si — s — ry/A(s, si, g 2 



(16) 



with g?i )2 = q 2 — 2pi^q and c?i + rf 2 = Si + — s. The integration is performed in the center 
of mass of the initial particles. 

The differential cross sections for heavy scalar and pseudoscalar states have been cal- 
culated for the reaction p + p —>■ i + + i~ + n^ir, as a function of the s and of the total 
momentum carried by the pions. Since the coupling constant g\ and the partial width Tp p 
are not known, the results for the differential cross section are rescaled by these quantities. 
For the numerical application, we focus on the scalar resonance /o(980) (Figs [3]) and on the 
pseudoscalar resonance, with 7/(1405) (Fig. Hj). The electromagnetic nucleon form factors 
have been parametrized according to 8]. 

The two dimensional plots show a structured landscape, limited by the kinematical cut 
due to the physical accessible region. Peaks appear in the projection to the q 2 axis (inte- 
gration over s\), and are due to the vector meson resonances coming from enhanced photon 
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FIG. 3: (Color online) Double differential cross section for the reaction p + p — > M* + n^ir, 
M* = /o(980), as a function of q 2 and si. 



vector meson coupling included in the FFs model. The projection to the si axis (integration 
over q 2 ) shows a bump which is due to the excitation of the pp resonant state. 

The absolute value of the cross section is in general about two orders of magnitude smaller 
than in the case of one pion production, as calculated in Ref . [l( , but still measurable in the 
resonance region. 

To avoid double counting, we do not consider direct emission of a virtual photon from 
the resonance, as well as from charged component from its decay, as we imply that it is 
implicitly taken into account in the resonance width. 

Let us discuss now the possibility for the creation of an atomic state of pN type. In the 
case of a charged atom (N = n) the situation is similar to the deuteron case, with a weakly 
bound state. One may expect several narrow states with masses < 2M and binding energy 
of ~ 2 — 3 MeV. The experimental observation of such states has been reported {4]. 

The coupling constant gi related to the creation of the bound state (pn) can be related 
(following 7]) to the wave function of this state: 

TO 



gx = M" 3/2 |^(0)|, J \ip\ 2 d 3 r = 1, #(0) = -M J V(r)r^{r)dr, (17) 
where \l/(r) and r are the wave function and the radius of the pn bound state and V(r) is 
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FIG. 4: (Color online) Double differential cross section for the reaction p + p — > M* + n^ir, 
M* = 17(1405) as a function of q 2 and s±. 

the potential of the p — n interaction. 

The reaction p + n — > 7Ti + ... + 7r n7r can be calculated using a similar formalism with 
evident modifications, replacing the proton FFs with the neutron FFs and adding a contact 
term proportional to (F™ — Ff)q ll /q 2 . One may expect that the general behavior of the 
distributions will be similar to the case of pp collisions. The size of the relevant resonance 
should be be much smaller than for pp one (several fermi), as for a deuteron-like object. 

Neutral atoms, (N = p) are of electromagnetic nature, their typical size can be estimated 
to 25 fm, and their binding energy of the order of 20 keV. The width of these resonances 
in S state is rather large. Such resonances, with orbital momentum L = have typically 
strong interaction width of order of 100 MeV. For 1^0, the width should be smaller, see 
Eq. (HJ). 

In conclusions, we have calculated the double differential cross section for the multipion 
production in proton antiproton collision, with emission of a leptonic pair. Numerical esti- 
mations show that the cross section is measurable in the kinematical region which will be 
accessible at FAIR and will allow to investigate the formation and the resonant structure of 



S 



a possible NN system. 
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